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We compute charmonium spectral functions in 2-flavour QCD using the maximum 
Q\ ! entropy method and anisotropic lattices. We find that the S-waves (J/ip and r] c ) survive 
up to temperatures close to 2T C , while the P-waves (x c o an d Xci) me h away below 1.3T C . 

OO : 

O . 1. INTRODUCTION 

■ 

The fate of charmonium at high temperature has generated much interest the past 
20 years [[!]. Contrary to the original prediction that J/ip production would be highly 
suppressed immediately above the phase transition, recent theoretical and experimental 
results indicate a more complicated picture, where the IS states may survive up to high 
temperatures (possibly even with enhanced yields), while other states disappear earlier. 

The properties of hadrons at high temperature are encoded in the spectral function pr, 
^ . which can be related to euclidean-time correlation functions G? as follows, 



, du . . coshL>(r — 1/2T)] 
J 2ir sinh(u;/2T) 

where the subscript T correspond to the different quantum numbers. Determining the 
continuous function p(uj) from G(r), r — 1 . . . L T , is an ill-posed problem, but it is possible 
to infer the most likely p(u) using the maximum entropy method (MEM) [ 2]. This has 
been used to study spectral functions in the quenched approximation [ El IH E], but 
substantial uncertainties still remain. 

In order to resolve some of these uncertainties, lattice simulations with dynamical 
fermions (2 or 2+1 light flavours) would be highly desirable. However, in order to deter- 
mine spectral functions using MEM, at least O(10) independent lattice points are needed 
in the imaginary time direction. At T ~ 2T C , this implies a temporal lattice spacing 
a t < 0.025fm. If the spatial lattice spacing a s were to be the same, such a simulation 
would be far too expensive to carry out with current computing resources. 

In order to make the simulation feasible, anisotropic lattices, with a t <C a s , are therefore 
required. However, dynamical anisotropic lattice simulations introduce additional com- 
plications not present in isotropic or quenched anisotropic simulations. The anisotropic 
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formulation introduces two additional parameters (the bare quark and gluon anisotropies) , 
which must be tuned so that the physical anisotropies are the same for gauge and fermion 
fields. In the presence of dynamical fermions, this requires a simultaneous two-dimensional 
tuning, which has been described and carried out in Ref. [Ej. 

We have previously presented preliminary results [EHB], indicating that S-wave states 
survive well into the plasma phase. These simulations were carried out before the tuning 
was completed, and therefore had systematic uncertainties at the 20% level. Here we 
present results from a simulation where the anisotropies have been tuned to within 3%, 
and study P-waves in addition to S-waves. 

2. RESULTS 

We use the TSI gauge action and a coarse Wilson, fine Hamber-Wu fermion action 
[ [TU] with stout-link smearing [ |TTJ as described in [ E] . The spatial lattice spacing is 
a s pa 0.18fm and the anisotropy £ = 6, while m n /m p = 0.54; our simulation parameters 
correspond to run 6 in [E]- For the charm quarks we have used bare masses a^m = 0.08 
and 0.092. The critical temperature was determined from studying the variation of the 
Polyakov loop as a function of N T on 12 3 x iV T lattices, giving N£ = l/aT c pa 34. 500 
configurations were generated on 8 3 x 32, 8 3 x 24 and 8 3 x 16 lattices, corresponding to 
T/T c pa 1.05, 1.4 and 2.1 respectively. Charmonium correlators were computed using all- 
to-all propagators [ 12J and analysed with Bryan's MEM algorithm [ using the free 
continuum spectral function u 2 as default model. 

Figure show spectral functions for the pseudoscalar and vector channels, correspond- 
ing to the S-wave states r\ c and J / ip respectively. At the two lower temperatures (N T = 32) 
and N T = 24) we see two peaks; however, the second (higher) peak is a lattice artefact 
which can be observed as a cusp in the free lattice spectral function [Ej. The lower peak 
position corresponds to the S-wave ground state masses at zero temperature, indicating 
that these states survive more or less unchanged at least up to 1.4T C . At T ~ 2.1T C the 
picture is less clear: for the lighter mass (m = 0.08) the peak has disappeared in both 
channels; while for the heavier mass (m = 0.092) there is still a peak in the vector channel, 
but not in the pseudoscalar. To what extent these effects are genuine, i.e. whether they 
show a real mass dependence and the survival of the vector meson to higher temperatures 
than the pseudoscalar, cannot be ascertained at present. 

In fig. |21 we show spectral functions for the scalar and axial- vector channels, correspond- 
ing to the P-wave states Xco and Xci respectively. In the axial-vector channel there is a 
peak at atOJ ~ 0.45 for the lowest temperature, consistent with the presence of a modified 
Xci at 1.05T C . There is no discernible corresponding peak in the scalar channel. At the 
higher temperatures the peak has disappeared, indicating that these states have melted 
below 1.4T C . The peak in the scalar spectral function at a t uo pa 0.25 at N T = 24 is not 
clearly understood at present and requires further investigation. 

3. DISCUSSION 

The main conclusion that can be drawn from this study is that the S-wave states J/ijj 
and T) c survive virtually unchanged in the medium above T c , before melting at T < 2T C . 
The P-wave states, on the other hand, disappear shortly above the phase transition. 




Figure 1. Spectral functions for the S-wave charmonium states r\ c (top) and J /ip (bottom), 
for a t m = 0.080 (left) and 0.092 (right). 



Our results thus lend support to the qualitative picture that has emerged from quenched 
calculations. In order to make these results quantitative, a number of outstanding issues 
needs to be addressed: 

• The 8 3 lattice corresponds to a physical volume of (1.4fm) 3 , and finite volume ef- 
fects are expected to be substantial, especially for the P-waves. We are currently 
repeating the calculation on 12 3 lattices. 

• A more detailed temperature scan will be necessary to determine the melting points 
of the different states. This is currently underway. A zero-temperature run is also 
in progress to provide a baseline for comparison. 

• The present analysis has been carried out using only the free continuum spectral 
function as default model. In order to assess systematic uncertainties arising from 
the MEM, different default models will need to be employed. In particular, the free 
lattice spectral function [|7j will provide a useful comparison. 

In the final instance, simulations on finer lattices will be necessary to resolve the remain- 
ing uncertainties. This will however require a new nonperturbative tuning exercise, and 




Figure 2. Spectral functions for the P-wave charmonium states Xco and Xci, for «t^o — 
0.080 (left) and 0.092 (right). Thick lines correspond to Xco, thin lines to x&- 



remains a longer-term prospect. 
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